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Abstract 

Three bosons with large scattering length show universal properties that do not depend on the 
details of the interaction at short distances. In the three-boson system, these properties include 
a geometric spectrum of shallow three-body states called "Efimov states" and log-periodic depen- 
dence of scattering observables on the scattering length. We investigate the modification of the 
Efimov states in a finite cubic box and calculate the dependence of their energies on the box size 
using effective field theory. We explicitly verify the renormalization of the effective field theory in 
the finite volume. 
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Strongly interacting quantum systems can show universal properties that are independent 
of the details of their interaction. These properties establish connections between systems 
over a wide range of scales and with different underlying interactions. One example are the 
hydrodynamic properties of quantum liquids. The strongly interacting quark gluon plasma 
created at the Relativistic Heavy Ion Collider and ultracold Fermi gases both behave as 
nearly perfect liquids with almost no viscosity [l| . Their viscosity to entropy density ratio 
is close to a lower bound that was conjectured using string theory methods [2J. 

Another example are universal properties in non- relativistic few-body systems [jjj]. Here 
the strong interaction regime is characterized by a large scattering length a. If a is large 
and positive, two particles of mass m form a shallow dimer with energy E2 ~ —h 2 /(ma 2 ) , 
independent of the mechanism responsible for the large scattering length. Examples for such 
shallow dimer states are the deuteron in nuclear physics, the 4 He dimer in atomic physics, 
and possibly the new charmonium state X(3872) in particle physics 0]. These systems span 
more than 13 orders of magnitude in energy ranging from MeV to neV. In the three-body 
system, the universal properties include the Efimov effect [I| . If at least two of the three pairs 
of particles have a large scattering length \a\ compared to the range r of their interaction, 
there is a sequence of three-body bound states whose energies are spaced geometrically 
between —ff/mr 2 , and —h 2 /ma 2 . In the limit 1/a — > 0, there are infinitely many bound 
states with an accumulation point at the three-body scattering threshold. These Efimov 
states or trimers have a geometric spectrum Q]: 

E*f = -(e- 2 *l s °) n - n *h 2 Kl/m, (1) 

where is the binding momentum of the Efimov trimer labeled by n*. This spectrum is 
a signature of a discrete scaling symmetry with discrete scaling factor e n ^ s °. In the case 
of identical bosons, Sq ~ 1.00624 and the discrete scaling factor is e n ^ s ° ~ 22.7. This 
discrete scale invariance is also relevant if a is large but finite. It becomes manifest in 
the log-periodic dependence of scattering observables on the scattering length a Q. The 
consequences of discrete scale invariance and Efimov physics can be calculated in an effective 
field theory for short-range interactions, where the Efimov effect appears as a consequence 
of a renormalization group limit cycle [6j. 

Experimental evidence for an Efimov trimer in ultracold Cs atoms was recently provided 
by their signature in three-body recombination rates Q. This signature could be unrav- 
elled by varying the scattering length a over several orders of magnitude using a Feshbach 
resonance. More recently, evidence for Efimov trimers was also obtained in atom-dimer 
scattering [8[ and possibly in three-body recombination in a balanced mixture of atoms in 
three different hyperfine states of 6 Li [9], [Io| . 

The observation of Efimov physics in nuclear and particle physics systems is complicated 
by the inability to vary the scattering length. See, e.g., Ref. [llj for a recent discussion of 
Efimov physics in halo nuclei. Another opportunity to observe Efimov physics is given by lat- 
tice QCD simulations of three-nucleon systems [l2j. A number of studies of the quark- mass 
dependence of the chiral nucleon-nucleon (NN) interaction found that the inverse scattering 
lengths in the relevant 3 Si- 3 D! and 1 Sq channels may both vanish if one extrapolates away 



from the physical values to slightly larger quark masses [13J, [14j, |15j. This implies that QCD 
is close to the critical trajectory for an infrared RG limit cycle in the three-nucleon sector. 
It was conjectured that QCD could be tuned to lie precisely on the critical trajectory by 
tuning the up and down quark masses separately [161]. As a consequence, the triton would 
display the Efimov effect. More refined studies of the signature of Efimov physics in this 
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case followed [171 . |18| |. However, a proof of this conjecture can only be given by an obser- 
vation of this effect in a lattice QCD simulation [l2|. The first full lattice QCD calculation 
of nucleon-nucleon scattering was recently reported in 19[ but statistical noise presents a 



serious challenge and no three-nucleon calculation has been carried out to date. Since lattice 
simulations are carried out in a cubic box, it is important to understand the properties of 
Efimov states in the box. Apart from this application, the modification of Efimov physics 
in a finite volume is an interesting question on its own and could be tested with ultracold 
atoms in optical lattices. 

The modifications to the Efimov spectrum can be calculated in effective field theory (EFT) 
since the finite volume modifies the infrared properties of the system. The properties of three- 
body systems in a finite volume have been studied previously. For repulsive and weakly 
attractive interactions without bound states, Tan has determined the volume dependence 



of the ground state energy of three bosons up to 0((a/L) 7 ) [20] . In Refs. [2l|, [22], this 



result was extended for general systems of N bosons. First results for three and more boson 



systems from lattice QCD have become available recently [23J, |24], [25j. For the unitary limit 
of infinite scattering length some studies have been carried out as well. The properties of 
three spin-1/2 fermions in a box were investigated in (26[. In this system, the Efimov effect 
does not occur due to the Pauli principle. The spectrum of three bosons in a harmonic trap 
was calculated by Werner and Castin [27j ] . 

In this letter, we investigate the finite volume corrections for a three-boson system with 
large but finite scattering length. In particular, we derive an integral equation for the trimer 
energies in a system of three identical bosons in a cubic box and solve this equation for a 



few exemplary cases. We follow the strategy proposed in [28( and extend the the EFT for 
the three-boson system with short-range interactions of Ref. ((J to a cubic box. We also 
verify explicitly that the renormalization of the EFT is not affected by the box. Our study 
is carried out to leading order in the large scattering length. This corresponds to the zero 
range limit with r$ = 0. Our results will be applicable to physical states whose size is large 
compared to tq. Using units with H — m — 1 from now on, the effective Lagrangian can be 
written as 0] 

C = ^ (id t + -V 2 ] if> + ^dU - ^ (efty 2 + h.c.) - ^dt#V + . . . ; ( 2 ) 
\ 2 J 4 4 36 

where the dots indicate higher order terms. It involves the boson field ip and an auxiliary 
dimer field d. The coupling constants #2 and g% are matched to the scattering length and a 
three-body observable, respectively. 

The atom-dimer amplitude is determined by the inhomogeneous integral equation de- 
picted in Fig. [TJ where the single lines denote the boson propagator and the double lines 
denote the full dimer propagator. In a cubic box, the momenta are quantized and the 
loop integrals in the equation are replaced by discrete sums over quantized momenta 
k = (2n/L) n, n G Z 3 , where L is the side length of the box. The loop sums are regulated 
via a momentum cutoff A just as in the infinite volume case. All observables are indepen- 
dent of A. This is achieved by an appropriate dependence of the coupling constant g% on A 
0]. This dependence can be calculated analytically, but one three-body input is required 
to fix the unknown constant A* (see below). This renormalization procedure is performed 
in the infinite volume limit. The ultraviolet physics associated with the renormalization of 
the EFT should not interfere with the infrared physics brought about by putting the system 
into a finite volume, as long as the corresponding scales are well seperated. We will verify 
this expectation explicitly below. 
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FIG. 1: Integral equation for the atom-dimer amplitude. The single lines denote the boson propa- 
gator and the double lines denote the full dimer propagator. 

The full dimer propagator D denoted by the double lines in Fig. [1] is obtained by dressing 
the bare dimer propagator which is simply a constant, 4z/g 2 , with bosonic loops. This leads 
to an infinite sum which can be evaluated analytically yielding 

-l 

(3) 

for a dimer at rest. The term containing the box length L vanishes in the limit L — > oo and 
the expression reduces to the infinite volume result. Using the Feynman rules encoded in 
Eq. (j2J) and the full dimer propagator from above, we can translate the Feynman diagrams 
in Fig. [1] into an equation for the atom-dimer amplitude. It involves an integration over 
the loop energy and a sum over over the quantized loop momenta. The integration over 
the loop energy is performed using the residue theorem while the remaining sum over the 
quantized momenta is rewritten into an integral by virtue of Poisson's resummation formula: 
Sriez 3 ^ 3 (-2 — ^) — J2m£Z 3 exp(i27rm • z) , which is understood to be used under an integral. 

The resulting expression is simplified further by exploiting the behavior of the atom-dimer 
amplitude near a bound state. For energies close to a bound state energy, the amplitude has 
a simple pole, and the dependence on the ingoing and outgoing momenta seperates. Thus, 
we obtain a homogeneous integral equation for the bound state amplitude J 7 , namely 

r(p) = \ f d3 y E ^^(p, y)r E (y)F(&) , (4) 

where 

Ze(p, V) = (p 2 + P ■ V + V 2 -E)-' + ^, (5) 

The A-dependent three-body interaction in Eq. ([5]) is given by H(A) = — (73A 2 / (9(7|) = 
cos[so k^A/A*) + arctanso]/ cosfso ln(A/A*) — arctanso] where so ~ 1.00624 is a transcen- 
dental number 0. 

In the infinite volume case, only S-wave bound states are formed. In a cubic box, however, 
the extraction of the S-wave part of T is not straightforward due to the breakdown of the 
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spherical symmetry to a cubic symmetry. The infinitely many irreducible representations 
of the spherical symmetry are mapped onto the five irreducible representations of the cubic 
symmetry. To investigate the mixing of the different partial waves, the part of the amplitude 
with quantum number £ is projected out 



Fe(p)c m =- dyy 



nei a t',m' t",m" 



£' t" i 
m' m" 



(2£' + 1)(2£" + 1) // 
x \l 2iTl ie ''h''{Mn\y)Y t/m o(n n )Z { E '(p,y)r E (y)c e/ml F e (y) 



where we have used Wigner 3-j symbols and 



2£+l 



i n f P 2 + y 2 -E \ H(A) 
py V vy J A 



(6) 



(7) 



with Qn a Legendre function of the second kind. Moreover, j#i is a spherical Bessel function 
of order £". The £" sum runs over all partial waves while the £' sum runs over partial waves 

associated with the A\ representation of the cubic group, namely £' = 0, 4, 6, The size 

of the contribution of the respective partial wave to the cubic group harmonic of the A\ 



301] ). The first term in 



representation is given by the coefficients cg m (See, e.g., Refs. [2£ 
Eq. reproduces the equation for the corresponding partial wave of the infinite volume 
amplitude, while the second term yields corrections due to the quantization of the momenta 
and admixtures of other partial waves. 

For our study of the trimer states in the cubic box, this equation is now specialized to 
the £ = case: 



71 Jo V f-? 3 L\n\y ) 



(8) 



where the dots indicate admixtures from higher partial waves. Since the leading term in the 
expansion of the Bessel functions in Eq. ([B]) is l/(L\n\y) these contributions are suppressed by 
at least a/L. They will be small for volumes not too small compared to the size of the bound 
state. The lowest partial wave that is mixed with the S-wave is the £ = 4 wave. Moreover, 
contributions from higher partial waves will be suppressed kinematically for shallow states 
with small binding momentum. This is ensured by the spherical harmonic in the second term 
of Eq. ([6]). Only for small lattices, i. e. when a/L is large, this behavior is counteracted 
by terms stemming from the spherical Bessel function jt(L\n\y) and higher partial waves 
may contribute significantly. We will therefore only consider the terms explicitly given 
in Eq. (jSJ) and leave the calculation of corrections due to higher partial waves for future 
work. For the numerical solution of Eq. (JHJ), a set of basis functions is chosen to turn 
the integral equation into a homogeneous matrix equation. The basis functions chosen are 
shifted Legendre polynomials defined on the intervall [0, A] with a logarithmic dependence 
in the argument. This behavior is tailored to the known log-periodic behavior of the infinite 
volume bound state amplitude [6|. The set of basis functions is orthonormal with respect to 
a suitable chosen scalar product. The energy E is now varied such that the resulting matrix 
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FIG. 2: Variation of the trimer energy E3 with the length L of the cubic volume for the states 
II (left) and III (right). Plotted are three datasets for different values of the cutoff parameter A, 
together with the l/(Ao) bands. The point a/L = corresponds to the infinite volume limit. 



has the eigenvalue one. In order to reduce the numerical effort, a Taylor expansion of the 
integral kernel 

4% VMV) = Z%(p,y)T Er {y) + A zP(p,y)r E (y) (E - E?) + O ((E - E?)*) 

i (111/ E=E^ D 

(9) 

around the infinite volume energy Ef is performed. In the following, we consider only the 
linear term in this expansion. We expect this approximation to be valid for small energy 
shifts; its validity will be discussed in detail below. More details on the numerical method 
will be presented in a forthcoming publication. 

For this first study of the trimer energies in a cubic box, we take a > and choose four 
states that have different energies in the infinite volume including shallow as well as deep 
states: 

la: £3°° = -1.18907/a 2 , A*a = 5.66, 

lb: £3°° = -27.4427/a 2 , A*a = 5.66, 

II: £3°° = -5.04626/a 2 , A*a = 1.66, 

III: £3°° = -11.1322/a 2 , A*a = 3.66, 

where is the trimer energy in the infinite volume. Note that the states la and lb appear 
in the same physical system characterized by A*a = 5.66. 

For each of these states, the energy in finite volume has been calculated for various L. 
Our results are shown in Figs. [2] and [3J In order to check the consistency of our results, the 
calculation was carried out for several cutoff momenta A and proper renormalization was 
verified explicitly. For each cutoff the three-body interaction H(A) was adjusted to give the 
correct energy in infinite volume. The calculated values in Figs. [2] and [3] agree with each 
other very well within the depicted error bands of the size l/(Aa) as long as the finite volume 
energy shift is not larger than 15-20% of the infinite volume energy. These error bands give 
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FIG. 3: Variation of the trimer energy £"3 with L for the states la (upper curves) and lb (lower 
curves). Plotted are two (three) datasets for different values of A, together with the l/(Aa) bands. 



an estimate of the finite cutoff corrections and renormalized quantities have to agree within 
these errors. They do not correspond to the error from higher order corrections in the EFT. 

When the finite volume shifts become larger than 15-20% of the infinite volume energy, 
the calculated results for different cutoffs in Fig. [2] start to disagree and proper renormal- 
ization cannot be achieved anymore. For the shallower state II this occurs for a/L ~ 0.6 
while the deeper state III can be treated in the linear approximation up to a/L rs 0.85. This 
behavior of the energy shift for the shallow states indicates the breakdown of the Taylor 
expansion of the integral kernel. As expected, the linear approximation in Eq. is not 
applicable anymore if the finite volume shift becomes to large. Calculations that do not 
rely on this expansion require a higher numerical effort. Such calculations are currently 
performed and the results will be presented in a forthcoming publication. 

In Fig. [31 we show our results for the two states la and lb in the same system characterized 
by A*a = 5.66 for box sizes where the linear approximation in the kernel is justified. For 
the shallowest state la reliable results can not be obtained for box sizes L < 3a while for the 
deepest state lb no appreciable finite volume shift can be detected for L > a. Generally, the 
L dependence of the energy shows a similar behavior for all investigated states. For large 
box sizes, the finite volume energy agrees with the infinite volume value. If the box size is 
decreased, the binding increases. The finite volume corrections are most important for the 
shallowest states which are largest in size and feel the finite volume first. For example, the 
ratio of the energies of the states la and lb decreases from 23.08 in the infinite volume limit to 
19.5 for L = 3a. Note that this ratio differs from the discrete scaling factor exp(27r/so) ~ 515 
even in the infinite volume limit. This behavior is expected for shallow states close to the 
bound state threshold [3J]. The ratio exp(27r/s ) will be approached if deeper states are 
considered. 

In this letter, we have studied the Efimov spectrum of three bosons in a cubic box with 
periodic boundary conditions. The knowledge of these modifications is important in order 
to understand results from future 3-body lattice calculations. Using the framework of EFT, 
we derived a general equation for the trimer energies in the box. We solved this equation for 
box sizes as small as L = a in a linear approximation of the kernel for small finite volume 
shifts and studied the modifications of the spectrum as the system is squeezed into the box. 
Moreover, we investigated the breakdown of the linear approximation. We find that the 
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expanded kernel seems to be applicable as long as the finite volume shift in the energy is not 
larger than 15-20% of the infinite volume energy. Numerically more expensive calculations 
using the unexpanded equation are under way. As a next step, the role of the higher partial 
waves contributing to the bound state amplitude has to be investigated in detail. This might 
allow to push the calculations to even smaller volumes and the regime L/a < 1 similar to 



the two-body system [281 ] . Finally, our method should be extended to the three-nucleon 
system in order to be able to test the conjecture of Ref. [l6[ in the future. This requires also 
the inclusion of higher order corrections in the EFT and finite temperature effects as lattice 
calculations will inevitably be performed at a small, but non-zero, temperature. 

This research was supported by the DFG through SFB/TR 16 "Subnuclear structure of 
matter" and the BMBF under contract No. 06BN411. 



[1] A. Cho, Science 319, 1180 (2008). 

[2] P. Kovtun, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 94, 111601 (2005) 

|arXiv:hep-th/0405231~] . 
[3] E. Braaten and H.-W. Hammer, Phys. Rept. 428, 259 (2006) |arXiv:cond-mat /0410417] . 
[4] V. Efimov, Phy. Lett. 33B, 563 (1970). 
[5] V. Efimov, Sov. J. Nucl. Phys. 29, 546 (1979). 

[6] PF. Bedaque, H.-W. Hammer, and U. van Kolck, Phys. Rev. Lett. 82, 463 (1999) 

|arXiv:nucl-th/9809025| ; Nucl. Phys. A 646, 444 (1999) | arXiv:nucl-th/9811046| . 
[7] T. Kraemer et al., Nature 440, 315 (2006). 
[8] S. Knoop et al., larXiv:0807.3306l [cond-mat]. 

[9] T.B. Ottenstein et al., Phys. Rev. Lett. 101, 203202 (2008) [arXiv:0806.0587 [cond-mat]]. 
[10] J.H. Huckans et al., arXiv:0810.3288 [physics. atom-ph]. 

[11] D. L. Canham and H.-W. Hammer, Eur. Phys. J. A 37, 367 (2008) [arXiv:0807.3258l [nucl-th]] 

and references therein. 
[12] K.G. Wilson, Nucl. Phys. Proc. Suppl. 140, 3 (2005) [arXiv:hep-lat/041204 3|. 
[13] S.R. Beane, PF. Bedaque, M.J. Savage, and U. van Kolck, Nucl. Phys. A 700, 377 (2002) 

[arXiv:nucl-th/0104030| . 

[14] S.R. Beane and M.J. Savage, Nucl. Phys. A 717, 91 (2003) |arXiv:nucl-th/0208021| ; Nucl. 

Phys. A 713, 148 (2003) |arXiv:nucl-th/0206113| . 
[15] E. Epelbaum, U.-G. MeiBner, and W. Glockle, Nucl. Phys. A 714, 535 (2003) 

[arXiv:nucl-th/02070"89l . 

[16] E. Braaten and H.-W. Hammer, Phys. Rev. Lett. 91, 102002 (2003) [arXiv:nucl-th/0303038] . 
[17] E. Epelbaum, H.-W. Hammer, U.-G. MeiBner and A. Nogga, Eur. Phys. J. C 48, 169 (2006) 
|arXiv:hep-ph/0602225] . 

[18] H.-W. Hammer, D.R. Phillips and L. Platter, Eur. Phys. J. A 32, 335 (2007) [ arXiv:0704.3726l 
[nucl-th]]. 

[19] S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, Phys. Rev. Lett. 97, 012001 (2006) 

[arXiv:hep-lat70 602010| . 
[20] S. Tan, Phys. Rev. A 78, 013636 (2008) [arXiv: 0709 .25301 [cond-mat]]. 

[21] S. R. Beane, W. Detmold and M. J. Savage, Phys. Rev. D 76, 074507 (2007) [ arXiv: 0707. 16701 
[hep-lat]]. 

[22] W. Detmold and M. J. Savage, Phys. Rev. D 77, 057502 (2008) [arXiv:0801.0763l [hep-lat]]. 



8 



[23] S. R. Beane, W. Detmold, T. C. Luu, K. Orginos, M. J. Savage and A. Torok, Phys. Rev. 

Lett. 100 (2008) 082004 [arXiv:0710.1827l [hep-lat]]. 
[24] W. Detmold, M. J. Savage, A. Torok, S. R. Beane, T. C. Luu, K. Orginos and A. Parreno, 

Phys. Rev. D 78, 014507 (2008) [arXiv: 0803 .27281 [hep-lat]]. 
[25] W. Detmold, K. Orginos, M. J. Savage and A. Walker-Loud, Phys. Rev. D 78, 054514 (2008) 

|arXiv:0807. 18561 [hep-lat]] . 
[26] L. Pricoupenko and Y. Castin, J. Phys. A 40, 12863 (2007) [ arXiv: 0705 .15021 [cond-mat]]. 
[27] F. Werner and Y. Castin, Phys. Rev. Lett. 97, 150401 (2006) l |arXiv:cond-mat/0507399] ; Phys. 

Rev. A 74, 053604 (2006) [arXiv:cond-mat/060782l1 . 
[28] S. R. Beane, P. F. Bedaque, A. Parreno and M. J. Savage, Phys. Lett. B 585, 106 (2004) 

[arXiv:hep-lat/0312004] . 
[29] F. C. von der Lage and H. A. Bethe, Phys. Rev. 71, 612 (1947). 
[30] S. L. Altmann and A. P. Cracknell, Rev. Mod. Phys. 37, 19 (1965). 



9 



